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Context. To expand the range in the colour-magnitude diagram where asteroseismology can be applied, we organized a photometry campaign 
to find evidence for solar-like oscillations in giant stars in the globular cluster M4. 

Aims. The aim was to detect the comb-like p-mode structure characteristic for solar-like oscillations in the amplitude spectra. The two dozen 
main target stars are in the region of the bump stars and have luminosities in the range 50-140 Lq. 
CN| Methods. We collected 6160 CCD frames and light curves for about 14 000 stars were extracted. The frames consist of exposures in the Johnson 
B, V and R bands and were obtained at three different telescopes. Three different software packages were applied to obtain the lowest possible 
photometric noise level. The resulting light curves have been analysed for signatures of oscillations using a variety of methods. 
J — Results. We obtain high quality light curves for the K giants, but no clear oscillation signal is detected. This is a surprise as the noise levels 
achieved in the amplitude spectra should permit oscillations to be seen at the levels predicted by extrapolating from stars at lower luminosities. 
|L* ' In particular, when we search for the signature of oscillations in a large number of stars we might expect to see common features in the power 
. 5^ I spectra, but even here we fall short of having clear evidence of oscillations. 

Conclusions. High precision differential photometry is possible even in very crowded regions like the core of M4. Solar-like oscillations are 
5_j probably present in K giants, but the amplitudes are lower than classical scaling laws predict. The reasons may be that the lifetime of the modes 
are short or the driving mechanism is relatively inefficient in giant stars. 

Key words. Stars: oscillations - Stars: activity - Techniques: photometric ~ Methods: Observational - Stars: evolution 



1. Introduction 

Asteroseismology has made great progress in recent years due 
to improved observational techniques. Extremely stable spec- 
trographs have been designed with the aim to detect exo- 
planets through Doppler shifts. These instruments have allowed 
for the unambiguous detection of /?-modes in about a dozen 
solar-type stars (for a review see Kjeldsen & Bedding '2004). 
Although most results have been obtained for dwarf and sub- 
giant stars, some giant stars have been shown to oscillate: ^ Hya 
(G7III; Frandsen et al. 2002l l. e Oph (G9.5III; De Ridder et al. 

Send offprint requests to: S. Frandsen 

* Based on observations with the Danish 1.54 m on La Silla at the 
European Southern Observatory in Chile, the 1.5 m at Cerro Tololo 
Inter- American Observatory in Chile, and the SSO40" at Siding 
Spring Observatory in Australia. 



I2006al l. and t] Ser (KOIII; Barban et al. l2004b . Unambiguous 
detection of /9-modes has so far been restricted to field stars in 
our immediate neighbourhood. 

DoUinger et al. ( I2005I I found that field stars show radial 
velocity and photometric variability, which seems to increase 
with increasing luminosity. This could either be due to a com- 
bination of activity and granulation or pulsations, either self- 
excited (Mira-like) or driven by convective motions (solar- 
like). Similar evidence for increased variability has been found 
for stars in the globular cluster 47 Tucanae (Edmonds and 
Gilliland [T996] l. 

Stellar clusters present a highly interesting prospect, due to 
the large number of stars collected in a limited field of view 
and additional constraints from cluster properties like relative 
age, metallicity, and evolutionary stage. Gilliland et al. d 19931 1 
were the first to make a search for solar-like oscillations in the 
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main sequence F-type stars in the open cluster M67. However, 
based on their multi-site campaign using 4-m class telescopes 
they made no clear detection. Recently, another campaign was 
carried out on the same cluster (Stello et al. 2006b, 2007 ) , this 
time concentrating on stars on the lower part of the giant branch 
with expected larger oscillation amplitudes. 

Knowing that p-modes are excited in giant stars, we de- 
cided to study K giants in the globular cluster M4. The two 
main reasons are that M4 has a rich population of K giants 
(see Fig.|2]i and the amplitudes are expected to be substantially 
higher than for subgiants. The prospect is that if we can mea- 
sure the gross properties of the oscillations, like amplitudes and 
the large splitting characteristic for solar-like oscillations, we 
could potentially study an ensemble of giant stars in different 
evolutionary stages. 

According to the empirical calibration by Kjeldsen & 
Bedding j 19951 1 oscillation amplitudes scale as L/(M T^ff)- For 
the K giants in M4, taking the solar amplitude as 4.7 ppm 
(parts per million), we expected amplitudes in the range 400- 
1000 ppm. The noise level we expected to reach in the ampli- 
tude spectra for the bright giant stars was ^ 25 ppm at high 
frequencies. This would allow us to detect any excess power 
and possibly the comb-like structure of the p-modes in the am- 
plitude spectra for the brightest and least crowded targets. The 
mode frequencies unfortunately decrease with luminosity and 
approach a range where noise caused by extinction and trans- 
parency changes is difficult to eliminate in ground-based data. 



Table 1. Observing log for the M4 campaign in 2001. 
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images 
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Observatory 


B 


V 


R 


Start 


End 


SSO 
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June 10 
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June 1 
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Fig. 1. Distribution of nights for the three observing sites. The 
time baseline is 78 days. 

and R and a few minutes in B depending on seeing conditions. 
Exposure times at SSO were 2-4 minutes in V and 1-2 minutes 
in R. On average this led to a mean duty cycle around 50%, 
since the readout time was 1-2 minutes. 



2. The observations 

CCD frames were obtained with the Danish 1 .54 m telescope 
at La Silla, the 1.5 m telescope at Cerro Tololo Inter- American 
Observatory (CTIO), both in Chile, and the 1 m telescope at 
the Siding Spring Observatory (SSO) in Australia. The obser- 
vations took place over a period of almost three months in order 
to get the necessary frequency resolution. We started at La Silla 
on April 13, 2001 and got the last data on June 27, 2001. We 
obtained useful data on 48 nights out of an allocation of 63 
nights. Some nights had overlap between Australia and Chile 
which can be used to check the consistency of the photometric 
results. Table [Tlgives a list of the collected dataset. In total we 
have 6160 CCD frames of M4 (4014 from La Silla, 863 from 
CTIO, and 1283 from SSO). 

In order to identify the mode type, we decided to observe 
in three filters B, V and R. Amplitude ratios between different 
colours will indicate which type of modes are present. To en- 
hance the signal to noise ratio the data can still be combined 
using an approximate scaling (Eq. |8]l as phase differences of 
p-modes between filters are small. About equal time was al- 
located to the three filters giving more frames in R due to the 
shorter exposure times. 

The time distribution is represented in Fig. [1] which in- 
cludes 3 1 nights from La Silla, 6 nights from CTIO, and 1 1 
nights from SSO. 

The exposure times were adjusted to give the best result 
for stars in the range 12.5 < V < 13.5, which is the area of 
the RGB bump in M4 shown in Fig. |2] Typically this meant 
exposure times at La Silla and CTIO of less than a minute in V 



3. CCD photometry techniques 

The images were first calibrated by performing traditional bias 
subtraction and flat fielding. In order to avoid nightly offsets, 
calibration frames were constructed for periods of several days, 
typically a week, and the same calibration frames were used for 
all images in that period. This point is quite important since the 
timescale of variability for the giant targets is in the range 9-3 1 
hours (cf. Sect. 5.1). 

CCD non-linearity was considered, and in the case of the 
La Silla data the non-linearity was corrected by the technique 
described by Stello et al. (I2006bl l. The other cameras did not 
provide the same option for a correction and were assumed to 
be linear. We have not seen any evidence for non-linearity, and 
as we were able to find comparison stars with similar magni- 
tudes in all cases, the influence of non-linearity should be neg- 
ligible. 

In Fig. [3] we present the calibrated CCD image, which has 
been used as the reference image for the La Silla observations 
of M4. The two dozen K giant stars we selected for a detailed 
analysis are marked by boxes. 

In order to achieve the highest possible photometric pre- 
cision we decided to use three different reduction programs, 
which have all been successful in the past in producing high 
quality results: isis, daophot and momf. In the most crowded 
areas we expected isis (Alard & Lupton ,1998. Alard ,2000 J to 
give the best results. In the less crowded areas we expected that 
DAOPHOT (Stetson 1 19871 1 would do best. Finally, we know that 
MOMF (Kjeldsen & Frandsen 1 1992| l performs very well in fields 
with only mild crowding. AU frames were reduced with isis 
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Fig. 2. The colour-magnitude diagram for M4 showing the bump stars and part of the horizontal branch. The colours and 
magnitudes have been calibrated using the database provided by Stetson (2007 ). To show the range covered, some of the stars 
selected for a detailed analysis are indicated and labeled with the ID - see Table|2] 



and DAOPHOT and a subset of the V frames from La Silla were 
reduced with momf. Each of the three methods will be described 
in Sects. [3n43.3l and we compare the results in Sect. 13.41 

3.1. The image subtraction method fisisj 

The main advantage of the difference-image technique is the 
ability to extract high precision photometry in the crowded re- 
gions near the core of M4. In addition, any variations due to 
airmass and transparency changes are removed to first order as 
a part of the image subtraction. 

The image subtraction method isis (Alard & Lupton [T998l 
Alard 120001 ) was used by two members of the team indepen- 
dently (isisl and isis2). For each filter and each observing site 
we selected the images with the best seeing to make the ref- 
erence images. While isis2 only used a single reference image 
for each site and filter, isisl used subsets of data from about one 
week. 

For each observed image, isis computed a kernel which de- 
scribes the variations of the PSF across the image relative to the 
reference image, isis then convolves the reference image with 
the kernel and subtracts this from the observed image. The re- 



sulting difference image will contain the signal that is intrinsi- 
cally different from the reference image, e.g. cosmic ray hits, 
hot pixels, and variable stars. 

In the first approach (isisl) only La Silla and CTIO data 
were included (Bruntt |2003l ). They were divided into five sub- 
sets and reduced separately. The reason for this was large dif- 
ferences in the positional angle of the CCD camera between 
different periods. We had problems using the photometry pack- 
age that is part of isis. Instead, we modified the aperture pho- 
tometry routine in daophot in order to be able to run it on the 
subtracted images, in particular allowing negative flux values. 
The necessary modifications were described in detail by Bruntt 
et al. ( 120031) . isis needs a reference flux for each star, as only 
the change in the flux is calculated. This was taken from the 
DAOPHOT reduction of the reference images. 

In the second reduction approach (isis2) we constructed 
only one reference frame (for every observing site and pass- 
band). The original isis software does not work correctly with 
frames that are significantly rotated with respect to the refer- 
ence frame. The problem resides in the interpolation process. 
We used our own application, which performs this task cor- 
rectly even for large rotation angles but requires much more 



4 



S. Frandsen et al.: A search for solar-like oscillations in K giants in the globular cluster M4 



I I I 



I I I 



I I I 



I I I 



I I I 



■ 78-':\..; 



■■243' . 



u 
o 



T 1 1 1 1 1 1 1 1 1 — I 1 1 1 1 1 1 1 1 r 

-6 -4-2 2 

arcmin 



T 1 r 



6 



Fig. 3. The field of M4 observed with the Danish 1 .54 m telescope at La Silla. The two dozen K giants we analysed in detail are 
marked by squares. The center of the field is at aaooo - 16''23™33''.7, ^2000 - -26°31'05".7. 



computation time. For all stars detected in the reference frames 
(of all three observing sites) we derived dififerential fluxes from 
our CCD frames following standard procedure of isis reduc- 
tions (for details see Kopacki .2000) . We performed both the 
PSF fitting and aperture photometry on the difference frames. 
To transform differential fluxes into magnitudes we used total 
fluxes measured in the reference frames using aperture pho- 
tometry tasks NEDA-t-DAOGROw Under daophot. Finally, all the 
data were merged into a uniform magnitude system by apply- 
ing magnitude offsets determined from a carefully chosen set of 
bright unsaturated stars common to all three observing fields. 



In this way, the CTIO and SSO data were transformed into the 
magnitude scale of the La Silla data. 

3.2. daophot/allstar/allframe reduction 

For the photometry we used the suite of photometry programs 
developed by Stetson ( |1987||1990l[1994 ): daophot, allstar and 
ALLFRAME. The general use of these is well described in manu- 
als and the literature. For our applications we adopted a slightly 
modified procedure, which consisted in the following. Firstly, a 
few of the deepest and lowest-seeing frames were used to pro- 



S. Frandsen et al.: A search for solar-like oscillations in K giants in the globular cluster M4 



5 



0.100 ^ 




Fig. 4. The mis scatter around the mean for a good night with 
the Danish 1.54 m telescope for the isisl reduction. The (red) 
squares give the noise level for the final time series for the se- 
lected K giants, where we have looked for solar-like oscilla- 
tions. The line indicates the estimated noise from scintillation 
and photon statistics. The excess of stars with high noise near 

V = 13.3 + 0.3 are the RR Lyrae stars. Saturation sets in at 

V ~ 12.5. 
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Fig. 6. The mean rms scatter per data point per night is plot- 
ted against the nightly mean flux level in the exposures repre- 
sented by a corresponding magnitude difference between each 
night and a reference night. The filled symbols are two nights 
with very few observations. The over-plotted diamonds indi- 
cate nights with a high sky level and the triangles nights with 
very good seeing. The rms uncertainty is marked for each point. 
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Fig. 5. Same as Fig. ID but for a bad night. The night is charac- 
terized by very large sky background contribution due to moon 
light. 

duce a master list of stars, which was obtained by several itera- 
tions through DAOPHOT, ALLSTAR and ALLFRAME. Ncxt, all frames 
were run through daophot/allstar once. With the preliminary 
photometry in hand, positional transformations between each 
frame and our reference frame were derived using daomatch 
and DAOMASTER (kindly provided by P. Stetson). 

From the frames used to derive the master star list, we cre- 
ated two lists of stars for generating the PSF for the individual 
images. The first list contained ~10 stars that were well isolated 
and these were used to generate the first version of the PSF for 
each frame. Next the second list, containing ~200 stars cover- 
ing the observed field, was used to generate the final PSF. We 
used this PSF when running allframe on all the images to pro- 
duce the final photometry. We made a cross id of all stars in the 
photometry files and the 7,000 brightest stars were selected for 
further analysis. 



3.3. MOMF reduction 

As an independent test we used the momf package (Kjeldsen & 
Frandsen II 992l l on the La Silla images to derive time series in 
the V band. The reduction had to be done in a slightly different 
way than normal, as the La Silla images were slightly rotated 
relative to each other momf can only handle shifts. Thus we had 
to introduce an additional coordinate transformation changing 
the (x, y) coordinates for each image before doing photometry. 
For stars outside the very crowded regions, results were compa- 
rable to the other techniques. We conclude that we reach similar 
results for non-crowded regions by all of the techniques. 

3.4. Comparison of the photometry 

Our database contains photometry from the three different re- 
duction methods and comprises light curves of 13 611 stars 
with up to 6 160 data points. We also store information about 
each frame such as seeing, airmass, background level, the mid- 
time of each observation (Julian date), and position on the ref- 
erence frame of each star. The data can be accessed at the web 
page http://astro.phys.au.dk/~srf/M4/ 

The rms scatter in the light curves is plotted versus magni- 
tude for a relatively good and bad night in Fig.|4]and|5] respec- 
tively. These results are based on the isisl photometry in the 
V-band. For the bad night, the noise is considerably larger than 
the good night for all stars. The higher noise is due to high sky 
background as a result of the close proximity of the Moon to 
M4. The lower envelope of stars for the good night is close to 
the expected precision based on the flux level (indicated by the 
solid line). We find a number of stars with rms scatter as low as 
1 .5 mmag. For the K giants the level is typically about 2 mmag 
per data point. The bump stars are found around V - 13.5, 
where we obtain the lowest noise levels. 
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Fig. 7. The mean scatter per data point per night for night 14 for the daophot (left panel) and isisl reduction (right panel). The 
line gives the noise level for a good night, see Fig.|4] 



As shown in Fig. |5] some nights do not give good results, 
with noise levels considerably higher than the majority of the 
data. Naturally, the noise per data point depends on the level 
of exposure per frame. This will vary as we had to adjust the 
exposure time to avoid saturation. In good seeing we obtained 
fewer photons per frame as the stellar image size is smaller and 
the maximum exposure level is kept constant to avoid satura- 
tion. Also, when cirrus clouds were present we kept the ex- 
posures times a bit conservative in order to avoid occasional 
over-exposed images from holes in the clouds. 

We have investigated in which conditions the photometry 
was worse than expected from photon noise alone. To do this 
we considered the isisl photometry of stars from La Silla. We 
selected an ensemble of 19 bright stars and calculated their av- 
erage scatter per night. The average flux for the ensemble was 
also derived and converted to a magnitude. One good night was 
chosen as the reference night and the magnitude for that night 
was subtracted from the magnitudes from all other nights. 

The measured rms scatter is plotted against the magnitude 
for all the nights in Fig. |6l If only photon noise were present 
the scatter should follow the solid line in Fig. |6l ie. following 
the relation: 5 logjQ cr - m — m^gf. 

Most of the nights follow this law, but there are exceptions. 
Nights with a high sky background clearly stand out as less 
valuable for the purpose of achieving high photometric preci- 
sion (diamonds in Fig. |6l). A second problem arises, when the 
seeing is very good since under-sampling sets in (triangles in 
Fig.©. 

An important issue is whether any of the three reduction 
techniques are superior In general we find that the results 
are quite similar for the image subtraction, PSF fitting, and 
aperture photometry techniques. However, momf is clearly not 
suited for crowded fields and only for a small fraction of stars 
does it perform as well as isis or daophot. The quality of the 
DAOPHOT photometry is comparable to the isis results, although 
isis is superior in the most crowded areas. We also find the best 
results with isis for faint stars, especially on nights with bad 
seeing. This is shown in Fig.|2]where we compare the rms scat- 
ter in data from one bad night from La Silla using daophot (left 



panel) and isis (right panel). The rms noise of the K giant stars 
is 40% lower for the isis reduction. 

It is possible that with a more careful choice of PSF stars 
and other parameters involved in the daophot reductions, it may 
be possible to improve these results. The optimum choice of 
parameters for good observing conditions might not be the best 
for frames obtained under non-photometric conditions. 

A final illustration of the differences seen between the dif- 
ferent algorithms is presented in Fig. [8] Here the rms scatter is 
presented for three different reductions: isisl, isis2 and daophot. 
Each point is for a particular star and the symbols define which 
reduction technique was being used. It is evident, that none of 
the techniques produce the best results in all cases. It seems that 
the image subtraction method isisl has been most successful on 
average. 

The time series from each observatory show different char- 
acteristics due to different observing conditions and different 
telescope and detector setups. At the SSO the seeing was con- 
siderably worse leading to more extended stellar images. This 
leads to increased noise as a result of the higher influence of the 
background level, which is worse for fainter stars. This effect 
is shown in Fig. |9] The CTIO data are heavily influenced by 
the presence of the Moon bing close to the target field, particu- 
larly in the middle of the run, and this makes the data of rather 
poor quality. The results can be characterized by the rms noise 
per data point for a non-crowded K giants for the three sites on 
a good night: 1.6mmag for La Silla, 3.0mmag for CTIO and 
2.1mmagforSSO. 



4. Light curve preparation 

The outcome of the previous section is either absolute or dif- 
ferential magnitudes. In both cases it is necessary to make dif- 
ferential photometry by comparing each target star with a set 
of reference stars. The effects of extinction, instrumental drifts 
and other sources of non-stellar noise need to be minimized 
for us to be able to study and stellar signal at frequencies be- 
low 30yL(Hz corresponding to periods longer than 9 hours. The 
approach is described in the following Sections. 
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Fig. 8. The rms noise level in K giant stars for La Silla data from two different nights is shown in the left and right panel. We 
compare the noise level from the three reduction methods. It is similar for most stars, but in a few cases the daophot and isis2 
noise levels are higher than for isisl. 



4.1. Ensemble photometry 

An ensemble of reference stars is used to derive the relative 
magnitude for all stars. This should take out transparency and 
extinction variations. In M4 there is a wide range of possibil- 
ities for choosing reference stars. However, the results depend 
strongly on the implementation of the ensemble average per- 
formed. 

The algorithm we employed was the following: for a given 
star, a number of reference stars were selected which were re- 
quired to have similar V magnitude and B - V colour and be 
within a certain distance (eg. half the size of the CCD) on the 
CCD frame. We calculated the magnitude difference of the tar- 
get star and each reference star for a whole observing block 
(t > 7d) and subtracted the median value producing a time 
series Aini{t), where the index / indicate the reference star. 

We now do an analysis night by night. For each night the 
rms noise o", of each Am,(f) time series is calculated and a 
weight is assigned using 
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The 10 stars with highest weights are used to calculate the final 
light curve for the night as the weighted sum 

10 10 

Aw(0 = Yj / 2 ^' • 

i=l 1=1 



4.2. Decorrelation 

Instrumental drifts and reduction noise may lead to correlation 
of the measured magnitudes with parameters like the colour of 
the star, position on the CCD and airmass. We found no clear 
correlations between the relative magnitudes and parameters 
Uke airmass, sky background or position on the CCD. As a con- 
sequence no type of decorrelation was performed. 



The set of reference stars will vary from night to night. The 
method does not imply any high-pass filtering at a timescale of 
one day, as the median is subtracted before splitting the time 
series into single nights. The technique gives slightly lower 
(10-20%) white noise levels in the final light curves compared 
with the simpler choice of a common set of reference stars ex- 
tending over all nights of an observing period. The final light 
curves are obtained by selecting, for each observing site, the 
result with lowest noise among the different methods applied 
(daophot, isisl and isis2). 



4.3. Combining B,V and R 

The results from the procedures described in Sect. 14. H are a set 
of light curves from the three sites in three filters B, V and R. In 
our search for solar-like oscillations in the K giants we assume 
that phase shifts between observations in different filters are 
small (Jimenez et al. |1999l l. We also assume that the oscillation 
amplitudes scale with inverse wavelength as A oc A^^ (Eq. [8j 
and we therefore scaled the B and R time series relative to the 
V filter: scaling factors were 1.222 and 0.846. The weights, de- 
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Table 2. The 24 K giants selected for detailed analysis. ID2 is from Stetson ( |2007| ) if present in the catalog (NA means not 
available) and V and B - V are calibrated using this database. T^ff is found using the B - V calibration of Alonso et al. dl 9991 1 
while the luminosity and radius were determined using the parameters for the M4 cluster (see text for details). In the last four 
columns we list the predicted values for the peak amplitude in ppm when assuming 6L/L oc (L/M)"^, v^ax is the expected 
peak frequency for ;?-modes, Avq is the large separation, and cjc is the acoustic cutoff frequency (frequencies are all in //Hz and 
calculated from Eqs.|5]-[8]l. 



ID 


ID2 


Q'2000 


<52000 


V 


B 


- V 








6L/L 




Avo 




66 


364 


16 23 53.65 
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scribed in the next subsection, were also scaled with the inverse 
factor following Eq.lH 

4.4. Removing outliers and weighting ttie data 

Due to the changes in observing conditions there is a large vari- 
ation in the data quality from night to night. It is therefore es- 
sential to apply weights when the Fourier spectra are calcu- 
lated. 

In the first step we remove obvious outliers. All points that 
deviate more than 5cr from the mean are discarded. Weights are 
then calculated for each star by combining a weight for each 
night and a weight based on the point-to-point scatter. Each 
night gets a weight, which is 

= l/o-^g^P (3) 

where cr„,g/,f is the rms noise of the time series for the entire 
night. AH data points from this night get the same weight. 

The point-to-point weight is calculated by first computing a 
spline fit to an averaged time series with a smoothing width of 
+0.075 d. The spline fit is subtracted from the data. The noise 
level at any given point in time cTp, is then the average rms 
calculated using a Gaussian with a width of 0.05 d. The weight 
is chosen as the inverse of the noise level: 

W2 = l/o-p,. (4) 



We note that Eq.|4](rather than use of inverse variance) was 
also the preferred weighting in similar multi-site campaigns by 
Handler (12003]) and Bruntt et al. (l2007l) . 

The final weight we use is the product of wi and W2- A 
maximum weight (4 times the mean weight) is defined to avoid 
nights with very few points to give very high weights. 

5. Selection of K giant targets 

We selected two dozen K giants for a detailed analysis. The 
stars are found in the least crowded parts of M4 and cover a 
range in luminosity on both sides of the bump stars. Further, 
we selected the stars with the rms scatter in the time series less 
than 4 mmag. We list their ID numbers, standard V magnitude 
and B - V colour in Table |2] The position of some of the stars 
in the colour-magnitude diagram is shown in Fig. |2] and the 
position in the cluster are marked with squares in Fig. [3] 

To calculate L/Lq, T^ff and R/Rq for the targets we used 
cluster parameters from Ivans et al. ( II999I I for M4: distance 
d -2.1 kpc, reddening E{B -V) - 0.37, interstellar absorption 
Av - 4.0E(B - V), metal content [Fe/H]= -1.2, bolometric 
correction BC - -0.48 and mass M - 0.85 Mq. As pointed 
out by Ivans et al. (1 19991 1. the reddening is varying over the 
field of the cluster These parameters lead to a distance modulus 
of (m - M)v = 13.1. This compares reasonably well with the 
mean magnitude of the RR Lyrae stars (Kopacki & Frandsen, 
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I2007I I in the range 13.2-13.3, which have absolute magnitude 
My - 0.0. The uncertainty of the cluster parameters introduce 
a ~1Q% uncertainty in the luminosities. 

5.1. Expected signal for solar-like oscillations 

We have made estimates of the amplitude and frequency that is 
expected for solar-like oscillations. 

We use the following scaling relations for the frequency 
separation Avo, the expected location of p-mode maximum fre- 
quency Vmax and the acoustic cutoff frequency a>c: 



Avo = 134.9(M/M0)'/2(/?/Ro)"^''^juHz, 



3050 



OJr = 5600 ■ 



(/?/Ro)2(reff/5777K)0-5 
(M/Mo) 



jt/Hz, 



jt/Hz. 



(5) 
(6) 

(7) 



The expected amplitude of the oscillation modes were esti- 
mated from the calibration by Kjeldsen and Bedding ( I1995I I. 
but using the scaling 6L/L oc (L/M)^ '' suggested by Samadi et 
al. (12006b : 



0.7 



(6L/Lh = 4.7 



[(L/Lo)/(M/Mo)] 
(^/550nm)(reff/5777K)2 



ppm. 



(8) 



The computed parameters using Eq. |5]-[8] for the 24 selected 
K giants are listed in Table |2] As the luminosity decreases 
from about 150 to 50 L© the amplitudes decreases from 300 to 
1 00 ppm, while the frequency of maximum power shifts from 9 
to 30;uHz (periods of 31 to 9 hours). The search for the signa- 
tures of the oscillations will be difficult, as we are in the region 
of the frequency spectrum that is affected by slow drifts origi- 
nating from observing conditions (seeing, transparency and ex- 
tinction) and the instrument (eg. slight changes in the position 
of the stars on the CCD). 



6. Time series analysis 

6.1. Fourier analysis 

The light curves were analysed using the Fourier analysis pro- 
gram period98 (Sperl 19981 1. In Fig. [TO] we show the power 
spectra computed from the light curves using weights (cf. 
Sect 14.4b . The spectra of a dozen K giants covering a range 
in luminosity increasing from 50 to 140 from the bottom to 
the top panels. 

In general the white noise increases as the stars become 
fainter but some stars have significantly higher noise than stars 
with similar brightness, which in some cases is due to the pres- 
ence of a close neighbouring star but in other cases may be 
an unidentified instrumental problem. The noise level in the 
amplitude spectra lies in the range 60- 100 ppm at 81-104/iHz 
(white noise) increasing to 100-140ppm at 23^6 yuHz. With 
this low noise level we might be able to detect oscillations if 
the amplitudes follow the extrapolations used in Sect. 15.11 (cf. 
Table |2]i and the mode lifetime is not as short as a few days 
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Fig. 11. The autocorrelation of the spectral window for star 
#181. The vertical lines indicate peaks at 11.57, 5.84 and 
0.81 /iHz (1, 0.5 and 0.07 c/d). Note that multiples of the 
0.81 fjHz peak are present. 



as found by Stello et al. (I2006al l. The stochastic excitation of- 
ten gives rise to higher than average amplitudes for observing 
periods that do not exceed the lifetimes by large factors. 

The observed amplitudes of the highest peaks below the 
acoustic cutoff frequency a>c range from 200 to 400 ppm in the 
spectra plotted in Fig.[TOl We cannot claim that they are real as 
the signal-to-noise is low (S /N < 4) in all cases. In a few stars 
we see suggestive evidence of power in the right range (#181, 
#228 and #229), which occurs between multiples of 1 c/d, but 
due to the low S/N we cannot claim to have detected p-modes. 

6.2. Quantifying the effect of the spectral window 

The observing window will introduce frequency spacings 
which may be mistaken for the frequency spacing characteris- 
tic for solar oscillations. The spectral window is different from 
star to star as the applied weights vary from star to star 

In Fig.[TT|we show the autocorrelation function of the win- 
dow function for star #181. This is a representative example 
for the K giant stars and gives an idea of the effect of the com- 
plicated spectral window which is due to gaps in the time se- 
ries. The autocorrelation function has a very pronounced peak 
at 11.57 yuHz due to the daily side lobes. In addition a strong 
peak at ~0.81//Hz is seen. The origin of the peak is found in 
the observing schedule, where dark time was not allocated to 
the project (Fig. [TJ. When searching for the presence of the 
large separation using the autocorrelation technique this will 
give rise to extra peaks. For a few stars the window function 
shows a more complex pattern giving rise to even more peaks 
in the autocorrelation. This will considerably reduce our ability 
to see clear evidence of stellar signal oscillations (cf. Sect. 17. lb . 



6.3. Subtracting 1 c/d aliases 

We note that some of the observed excess peaks in the am- 
plitude spectra are close to multiples of 1 c/d (11.57 fiHz) as 
indicated by the vertical dashed lines in Fig. \W\ These peaks 
are likely due to residual effects of extinction. In some light 
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Fig. 10. Power spectra for the 12 K giants with the lowest noise at high frequencies. The luminosity of the stars increases from 
the bottom to top panels and increases from left to right. The predicted location of the maximum p-mode power and the acoustic 
cutoff frequency (indicated by thick black lines) shift to lower frequencies for increasing luminosity. Frequencies at multiples of 
11.57yuHz (1 c/d) are indicated with vertical dashed hnes. 



curves a peak in the range v = 0.8-0.9 yuHz was also detected as 
likely caused by the spectral window (cf. Sect. 16.2b . Therefore, 
following the approach of Stello et al. ( I2006bl) . we subtract 
these low frequency peaks if they are significant {S /N > 4). 
We specify later in each case, whether we are using the raw or 
the cleaned spectra (cf. Sect.|7]i. 

In Fig. [I2]we show an example of the process for star #66. 
Daily alias variations are seen in the amplitude spectrum in 
panel (a). Panel (b) is the amplitude spectrum after subtracting 
two peaks indicated by the vertical tickmarks in panel (a). Panel 
(c) is the spectral window, which shows quite strong and com- 
plex sidelobes. The white noise level is ^ 60ppm and is calcu- 
lated as the mean level in the frequency range 81-104 yuHz. 



6.4. Light curve simulations 

To facilitate the interpretation of the observations in Sect. [T] we 
made simulations of the light curves. The time sampling and 



the noise properties are the same as for the observations and 
they include a known oscillation signal. 

To mimic the noise we used the observed time series of two 
stars assuming that they resembled pure noise. The two stars 
(#66 and #292) are at either end of the luminosity range of the 
selected K giants. Both are among the K giant stars with the 
lowest noise in the amplitude spectra. The two stars show no 
clear evidence of power excess or a frequency separation that 
could arise from stellar oscillations. 

We used different input amplitudes and life time of the 
modes using the simulation software described by De Ridder 
et al. (i2006b] l and SteUo et al. ( |2004t . 

7. A search for evidence of solar-like oscillations 

Inspection of the K giant amplitude spectra in Fig. [TOl reveals 
that the noise increases towards low frequencies, indicating 
that drift noise is present in the light curves. This will seri- 
ously hamper the interpretation of the signal below ^ 20 yuHz. 
We see no clear evidence of excess power or the comb-hke 
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Fig. 12. Amplitude spectra for star #66 which is the brightest, 
non- saturated star in the subsample chosen for a detailed anal- 
ysis. Panel (a) is the raw spectrum with one low frequency peak 
and a peak close to 2 c/d (23. 14 yuHz) marked. Panel (b) shows 
the amplitude spectrum after subtracting the two peaks marked 
with vertical lines in panel (a). Panel (c) shows the spectral win- 
dow, which is seen to contain a complex substructure around 
each 1 c/d sidelobe. 



pattern expected for solar-like oscillations in any of the K gi- 
ants. Although we find clear evidence for an increase in excess 
power towards low frequencies, it is not possible to disentangle 
drift noise and variation intrinsic to the K giant stars, eg. due to 
granulation or p-modes. 

We have used two techniques to look for evidence of p- 
modes in the amplitude spectra: 

- Autocorrelation of the amplitude spectra to search for the 
large separation (Sect. 17. it . 

- A search for evidence of excess power in the expected fre- 
quency range (Sect. |7.2l l. 

7.1. Autocorrelation 

A clear proof of the presence of stellar oscillations would be 
a systematic behaviour of the observed frequency separations. 
The separation should increase with decreasing luminosity (cf. 
Eq. |5]l. In a few cases the autocorrelation shows a peak that 
could be indicative of the large separation. A prominent peak 
is seen in star #181 as shown in Fig.[T3] where there is a peak 
at 2.75 juHz (predicted to be at 2.55 //Hz). The highest peak at 
^0.9 yuHz (vertical dashed line) is due to the window function 
(cf. Fig. [TTI ) The autocorrelation was done on the raw spec- 
trum. We find a few other stars with peaks in the autocorrela- 
tion spectrum that might come from the presence of a p-mode 
spectrum, but in the majority of cases we are not able to find a 
single prominent peak. 
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Fig. 13. The autocorrelation of the raw power spectrum for the 
giant star #181. Vertical lines represent peaks in the window 
function. The peak at 2.75 yuHz could be the large separation, 
but might also be due to low frequency noise interacting with 
the window function. 

We made several simulations of stars #66 and #292 with 
known input amplitude and lifetime as described in Sect. 16.41 
We used two methods to search for the known large separations 
of the inserted modes. The first method was a simple autocorre- 
lation technique and the second method cut up the spectrum in 
small pieces and added these pieces (like when an Echelle dia- 
gram is created) and then searches for peaks. The latter method 
is described by Christensen-Dalsgaard et al. ( 2007l l. Both meth- 
ods gave quite similar results. If the oscillations have a Rvalue 
above 200, the correct large frequency separation Avo is recov- 
ered in both stars for amplitudes down to 400 ppm. 

In Fig. [14] we show the autocorrelation of the spectrum of 
star #292 including simulated oscillations with large separation 
2.83yuHz, amplitude A =400 ppm and lifetime r=32d. We note 
that this lifetime is at least a factor two longer than expected. In 
the top panel the spectrum of #292 was cleaned for multiples 
of 1 c/d before adding the simulated signal. The autocorrelation 
clearly shows several peaks that correspond to the large separa- 
tion present in the simulated signal (red dotted lines). The first 
dotted line indicates the large frequency separation present in 
the simulation, which repeats several times. The blue, dot-dash, 
vertical lines indicate parasitic peaks in the window function at 
0.85 and 1 1.57/vHz (see Fig.fTTTi. In the bottom panel in Fig. [14] 
the procedure applied to the observed time series data was fol- 
lowed (see #181 in Fig.[T3]l. In this case the input large separa- 
tion is much harder to recover. 

For the bright star (#66) for inserted amplitudes as low as 
A = 200 ppm and large Q values (>30) the large separation 
is recovered. However, if the lifetime is shorter than ^16 days 
(Q < 12), we find that it is not possible to recover the input 
separation in any of the simulated series, even with amplitudes 
up to 800 ppm. Thus, our upper limits of 200^00 ppm only 
hold if the lifetime is greater than two weeks as predicted by 
Houdek & Gough (2002) and Houdek (l2006i) . and not a few 
days as found by Stello et al. (I2006al) . 



' Q is oscillation lifetime divided by period of the mode as used by 
Stello et al. ( I2006cb and presented for a set of stars in their Fig. 4 
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Fig. 14. Autocorrelation of a simulation of star #292 including 
oscillations. In the top panel the spectrum of #292 was cleaned 
for multiples of 1 c/d before adding the simulated signal and in 
the bottom panel, the same procedure was followed as for the 
observed time series data. The red dotted lines correspond to 
multiples of the large separation in the simulated signal. The 
blue vertical lines indicate parasitic peaks from the window 
function at 0.85 and 1 1.57yuHz (see Fig.fT2li. 



7.2. A search for evidence of excess power 

To search for excess power we followed the approach by Stello 
et al. (I2007I I used in their analysis of K giants in the open clus- 
ter M67. The K giant stars in M4 were put in three groups 
with luminosities around 60, 80, and 130 (faint, average or 
bright stars). For each group we calculated a mean raw ampli- 
tude spectrum. To demonstrate the effect of cleaning the spec- 
tra, we also produced a mean spectrum for the faint group of 
the cleaned spectra. 

In Fig. [15] we show the mean spectra. The raw spectra for 
the three groups of stars are quite similar with prominent peaks 
near multiples of 1 c/d. This is unexpected and suggests that the 
power seen possibly is dominated by a noise component due to 
instrumental drift in the data. 

Using the simulations described in Sect. l6.4l we probed our 
ability to identify the inserted oscillation signal. The black, 
solid spectrum in Fig. [15] is based on a mean of four time se- 
ries of star #292 which includes simulations of p-modes using 
different lifetimes in the range 2-8 days and peak amplitudes 
of A = 400 ppm per mode. This covers the range of expected 
lifetimes and averaging four simulations removes some of the 
stochastic features of the oscillations. This artificial spectrum 
is clearly above any of the observed mean spectra. Our con- 



300 
250 
200 
150 
100 

50 - 



1 



Faint, average and bright stars 
Simulated (A=400) L=75 




10 

Freq(/2Hz) 



100 



Fig. 15. Average smoothed spectra for three groups of ob- 
served K giants with luminosity L ^ 60,80 and 130 Lq. The 
faint group (raw and clean version) are the blue lines, the aver- 
age group is plotted as the magenta line and the bright group is 
the red line. The bottom curve is the faint group after multiples 
of 1 c/d are removed. The black solid line is for a simulation of 
star #292 with peak oscillation amplitudes of A = 400 ppm. 

elusion is that excess power from /?-modes with amplitudes 
A > 300 ppm would show up clearly in our observed spectra. 

8. Comparison with 3D-hydrodynamical 
simulations of granulation 

We have compared our results for the K giants with 3D hydro- 
dynamical simulations of granulation by Svensson & Ludwig 
(120051) . The effect of granulation is modelled by simulations of 
the convective motions in a box which is scaled to predict the 
effect expected to be observed when observing a spherical star. 
The result of Svensson & Ludwig (2005 ) indicate that K giants 
should have a low frequency component of granulation power 
similar to the noise levels achieved for some of the K giants in 
M4. 

One of their simulations (Ludwigl2006l kindly provided by 
G. Ludwig) is shown as the dotted black curve in Fig. [16] The 
stellar parameters are R/ Rq - 30.2, \ogg = 2.0, Tetf - 4560 K, 
and metallicity [M/H] - 0.0. This corresponds to a luminos- 
ity L = 385 Lq, which is larger by roughly a factor two than 
the brightest K giants we have observed. The larger luminosity 
means that the simulated granulation has larger amplitudes and 
excess power at lower frequencies than for our brightest target. 

A scaling of the solar power spectrum has been attempted 
by Kjeldsen & Bedding (in preparation) and by Stello et al. 
(I2007I I. A white noise component with the same level as the star 
chosen for comparison (#66) is added to a scaled granulation 
signal and a scaled p-mode spectrum. The scaling is very crude, 
and a difference of a factor two from our K giants would not be 
surprising. The granulation spectrum is shown in a logarithmic 
plot of power density in Fig. [16] 

The simulated spectrum (black dotted line) and the scaled 
spectrum (blue dashed line) both show a larger power density at 
low frequencies compared to the observations (red solid line). 
The peak of the observed spectrum coincides with the peak 



S. Frandsen et al.: A search for solar-like oscillations in K giants in the globular cluster M4 



13 



1 00.00 F 




0.1 1.0 10.0 100.0 

Freq(/iHz) 

Fig. 16. Power density for three cases: 1) the observed star #66 
(red, solid line), 2) a scaled power spectrum including granu- 
lation, p-modes with A - 300ppm and white noise as star #66 
(blue, dashed line), 3) the simulation by Svensson and Ludwig 
(l2005T l (black, dotted hne). 

around 10/iHz of the scaled spectrum, in which the power is 
due to the presence of oscillations. In the observed spectrum it 
is likely to be the 1 c/d alias. 

The reason for the high power at low frequencies in the 3D 
simulations is because the luminosity of the model simulated 
is twice the luminosity of star #66. The amplitude of the power 
density should scale with and the frequency with L ' , which 
will bring the 3D simulations closer to the other curves, but still 
above the observations. Svensson and Ludwig ( 2005) show that 
lowering the metal content in a 1 model makes the power 
density drop considerably. A change of metal content might 
remove the discrepancy between the observations and the 3D 
simulations of a giant. The 3D simulation (black dotted hne) 
falls below the other curves at high frequencies since no white 
noise was included. 

9. Conclusion and discussion 

We have presented results for a three-site multi-colour photo- 
metric campaign on the globular cluster M4. We used data in 
the Johnson BVR bands from 48 nights with a time baseline of 
78 days. 

We compared the precision of the photometry obtained 
from three different photometric packages. In the most crowded 
parts of the field we find the best results using a difference im- 
age analysis (Alard [2000l l. while classical aperture photometry 
give better results in the semi-crowded regions. The best pho- 
tometry is found for stars not affected by crowding, and in this 
case we achieve noise levels close to the theoretical noise limit. 

The main goal of the campaign was to detect solar-like 
oscillations in the rich population of K giant stars in M4. 
However, we are not able to claim an unambiguous detection 
and we summarize our conclusions here: 

- We detect no individual peaks with S /N > 4.0. 

- We find no significant large frequency separation from the 
autocorrelation of the amplitude spectrum. 



- We divided the stars with the best photometry in three 
groups and computed the average power spectra. From 
a comparison of the spectra with realistic simulations of 
granulation and an assumed comb-like distribution of p- 
modes, we can exclude that power is present from p- 
modes with peak amplitudes above 300 ppm whatever the 
damping time. This upper limit is in accordance with the 
(L/M)*^^ scaling suggested by Samadi et al. (I2006I I which 
predicts amplitudes below 300 ppm for the most luminous 
stars, but far below the scaling law L/M from Kjeldsen & 
Bedding ( fT995] l. 

- In all K giants the amplitude increases towards low frequen- 
cies. This would be consistent with granulation, but we find 
that the observed increase does not vary with stellar lumi- 
nosity. Thus we cannot differentiate between drift noise in 
the data and the expected granulation signal. 

- We have compared the observations with a 3D- 
hydrodynamical simulation for solar metallicity. The 
observed power density at low frequencies is much lower 
than the simulations of granulations predict. In their 3D 
simulations Svensson & Ludwig '2005' found that the gran- 
ulation power is substantially smaller for lower metallicity 
(M4 has [Fe/H] = -1.2). This could explain that our data, 
even for the brightest star, falls considerably below the 
simulations (Fig.[T6]l. 



Our negative results support the recent evidence for short 
lifetimes in giant stars (Stello et al. l2006al ) but contradict the 
theoretical predictions by Houdek & Gough (I2002I I. 

For the population II star v Indi, which has a low lumi- 
nosity (6 L0), p-modes have been detected (Bedding et al. 
I2OO6, Carrier et al. 2007). However, we cannot tell whether 
the stochastic driving still works in population II stars around 
L = 100 Lq, although we can say it cannot be very efficient 
and certainly not exceeding amplitudes following the (L/M)°^ 
scaling relation suggested by Samadi et al. (120061) . 

This is contrary to the conclusion by Stello et al. (2007} for 
the K giants in M67, who found indications of p-mode power in 
stars at lower luminosities (10-20 L0). They find better agree- 
ment between observations and predictions if the amplitudes 
scale as L/M. Their noise level at low frequencies do not per- 
mit any statement for stars at luminosities similar to the stars 
investigated in M4. It remains an open question, whether the 
stochastic driving still works without being destroyed by con- 
vective or radiative processes in population II K giants. 

Based on the results obtained from large multisite cam- 
paigns on M67 and M4, it seems difficult to get conclusive posi- 
tive results about the stochastic oscillations and the granulation 
from ground-based photometry. Improvements can be made, 
but extinction and other atmospheric effects will always be a 
limiting factor and especially so at the time scales of the K gi- 
ants. It is also difficult to get a clean window function due to 
changing weather patterns for campaigns lasting several weeks, 
which is unfortunately needed to study solar-like oscillations in 
K giant stars. 
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9.1. Directions for tlie future 

If a new campaign is organized we recommend to use a single 
filter, as colour information is not likely to be of a quality that 
can lead to mode identification. We found no use for the dif- 
ferences between the filters. We recommend that the V filter is 
used in the future: the B filter observations show strong colour- 
dependent extinction and the amplitudes will be very low in 
R. 

From the ground one would do better by organizing a 
longterm radial velocity campaign. One immediate advantage 
is that the noise from granulation is 10 times lower in velocity 
(Grundahl et al. 120071 their Fig. 1). One could observe several 
K giants during the night since the time scale of the oscillations 
is several hours. A dedicated network providing good time cov- 
erage would be needed in order to obtain a clean window func- 
tion. The proposed SONG network (Grundahl, et al. I2007I I is 
an example of such a programme. 

For fast rotating stars where the broad spectral lines prevent 
precise velocity results to be obtained, observations of solar- 
like oscillations can be made with photometry from space. 
This has resulted in the measurements of variability in pop- 
ulation I G giants by the WIRE satellite (Buzasi, Bedding & 
Retter l2003T l. and p modes in rj Boo with the MOST satellite 
(Guenther et al. lIOUSb . 

Giants are also among the targets for the new photometric 
satellite missions like the CoRoT mission (Baglin, Michel & 
Auvergne 120061) and the future Kepler mission (Christensen- 
Dalsgaard et al.l2007i The CoRoT mission observes the same 
part of the sky for up to 150 days while the Kepler mission will 
cover many more targets for up to six years. 
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